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(g) Attitude angle calculation apparatus for a geostationary satellite. 



(«) Attitude angle calculation apparatus is disclosed for 
calculating an attitude angle of a satellite very precisely. An 
earth sensor (3) detects the amount of deviation of an angle 
defined between a predetermined axis of the satellite and a 
straight line connecting the satellite and the geocentre of the 
earth. A fixed-star sensor (5) detects a physical intermatrix of 



a fixed star. A yaw angle determination device (4,6,7). in 
response to orbital position information (D) of the satellite 
sent from a ground station and outputs of the earth sensor 
and the fixed-star sensor, determines the yaw angle (iMt» of 
the satellite, independently of the roll angle (<b(t)| or the pitch 
angle (6(t)) thereof. 
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MITSUBISHI DENKI KABUSHIKI KAISHA 31/2176/02 
ATTITUDE ANGLE CALCULATION APPARATUS FOR A 
GEOSTATIONARY SATELLITE 

This invention relates to an attitude angle 
. calculation apparatus for use in precisely controlling 
the attitude of a geostationary satellite, more 
particularly an earth-directed satellite the attitude 
5 of which is controlled with respect to the earth in three 
orthogonal axes. 

For an earth surveillance satellite or a 
communications satellite, in order to maintain a 
predetermined orientation of an antenna or measuring 
10 device with respect to the earth, it is necessary to 

align a specified axis of the satellite with a straight 
line connecting the satellite and the earth. 

An attitude angle calculation apparatus is used to 
calculate the angular relationship of a coordinate system 
15 of the satellite with respect to a line connecting the 
satellite and the earth and a line orthogonal to a plane 
including the orbit of the satellite. An attitude control 
device controls the attitude angle of the satellite so 
that the result of the calculation performed by the 
20 attitude angle calculation apparatus becomes zero, to 
thereby maintain a predetermined axis of the satellite 
directed towards an object. 

Referring to Fig. 1 of the accompanying drawings, 
typically a conventional attitude angle calculation 
25 apparatus operates to detect relative deviation angles 
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$ and 6 of the Z g axis of the satellite / which is to 
be d i r ec t ed to the earth A, trojn the orbital 
coordinate system (Xq, Y^^^) measured by an earth 
whgor tn^untied op the pjttellite. In detail., the roll 
5 angle * r which is the angle of rotation of the Zg 
axis, of the satellite about its Xg axis from the Z Q 
axis of the orbital coordinate system, is detected 
by the earth sensor on the basis of the difference of a 
scanning width defined by fields of sight Bj and B 2 

10 of the earth sensor, and the pitch angle 9, which is 
thje. apgle 6f rotation of the Zg axis of the satellite 
about it* Y B axis» on the basis pf the amount of 
deviit^oij; of a reference pulse having a transmitting 
direction cy measured from the centre of the scanning 

15 ■ pulse fr However '. in- such £ conventional attitude 
angle <^W it is impossible to 

:.,/tnea$vrev:'i is the angle of . ' 

\tpt*ti6ii:)fjt$- the 0^ its- axis Zg.' 

: ^ it has been 

20 ^ U* by using a fixed-ptar 

|tett*oi^ Tig. 2 shows the concept 



- 3 - 



01 7471 5 



of such a proposal in which the star Polaris 2 (a-Urae 
Monalis) is detected by the fixed- star sensor directed 
along the -Yb axis of the satellite such ^that a center of 
the field of sight 1 of the fixed-star sensor is directed 
to Polaris 2. Assuming coordinates (Z,X) of Polaris 2 in 
the f ield of sight 1 of the fixed- star sensor, -the yaw 
angle y can be calculated as follows (provided that the 
satellite is in a stationary orbit) with the roll 
deviation and the pitch deviation always being zero, and 
the tilt of the orbit also being zero: 

0 = -(X + EpSin a) , ac[0, 2tr] , ... (1) 

where E p is the declination of Polaris 2 and a is the 
difference between the right ascension of Polaris 2 and 
the right ascension of the satellite. Since, however, it 
is difficult for a practical satellite to maintain the 
angles e and 0 and the tilting of the orbit at zero, the 
nonzero angles e and 0 and the orbital conditions of the 
satellite affect the measured value (Z,X) of the fixed 
star, and it is impossible to separate them from each 
other. Therefore, it is impossible . to measure the yaw 
angle ^ precisely. 

In another approach, a sun sensor is used as the 
fixed- star sensor. The sun sensor is mounted on the 
satellite with the field of sight thereof being in 
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coincid nc wij:h th -Xg axis of the satellite. In this 
case, the variation in of the angle of incidence n of the 
sun is detected. 

Fig. 3 shows the concept of the latter. in 
Fig. 3, reference characters 1,0, N, L, M, 8 and L* 
5 depict the orbital plane, the center of the earth A, a 
direction normal to the orbital plane, a cross line of a 
plane including a vector I in the direction of the sun and 
a line ON and the orbital plane 1, a point in the orbital 
plane 1 separated from the line L by 90*, an angle between 

10 a line OM and a line from the center 0 to a point P on the 
orbit, and a line parallel to the line L in a (X 0 ,Z 0 ) 
plane, respectively. 

As is clear from Fig. 3, when 6 ■ 0, i.e., the 
sun is always in the Xq,Yq plane, and assuming that the 

15 angle of incidence n is 0 at a certain time instant, the 
deviation a n thereof after the time instant coincides with 
a deviation &v of the yaw angle. That is, the roll angle 
fS is completely, separated from the deviation in . However, 
for a practical satellite, since the sun is out of the 

20 Xq,Yq plane due to the orbital movement of the satellite, 
the variation of the roll angle 4> about the Xq axis of 
the satellite affects in unseparably, causing precise 
measurement of yaw angle deviation *n to be impossible 
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(although there is no interference between the pitch 
angle 6 and An) . 

An object of the present invention is to provide an 
altitude angle calculation apparatus by which the yaw 
5 *angle can be obtained with high precision. 

According to the invention there is provided 
attitude angle calculation apparatus to be mounted on a 
satellite for calculating an attitude angle of the 
satellite, characterised by an earth sensor for detecting 

10 deviation of an angle defined between a predetermined 
axis of the satellite and a straight line connecting the 
satellite and the geocentre of the earth; a fixed-star 
sensor for detecting a physical intermatrix of a fixed 
star; and yaw angle determination means responsive to 

15 orbital position information of the satellite sent from 
a ground station and to outputs of the earth sensor and 
the fixed-star sensor to determine the yaw angle of 
the satellite independently of roll angle and/or pitch 
angle thereof. 

20 Embodiments of the invention will now be described, 

by way of example, with reference to the accompanying 
drawings, in which 

Fig, 1 illustrates the above-described concept of 
measurement of the roll angle and the pitch angle by an 
25 earth sensor; 

Fig. 2 depicts the above-described concept of yaw 
angle measurement by a fixed-star sensor; 

Fig. 3 shows a positional relation between the sun 
and a satellite and illustrates the above-described 
30 concept of interference between yaw angle deviation and 
the roll angle; 

Fig. 4 is a schematic block diagram of an embodi- 
ment of the present invention; 

Fig. 5 depicts a mathematical model showing the 



- 6 - 



0174715 



relationship between inertial space coordinates and 
orbital coordinates; 

Fig* 6 is a schematic block diagram of another 
embodiment of the present invention; 
5 Fig. 7 is a similar model to that shown in Fig. 5 

for use in explaining the embodiment of Fig. 6; 

Fig. 8 is a mathematical model showing the 
relationship between geocentric inertial coordinates and 
the position of the sun; and 
10 Fig. 9 is a schematic block diagram of a further 

embodiment of the present invention. 

Referring to Fig. 4 of the drawings, a first embodiment 
of an attitude angle calculation apparatus of the present 
invention includes an earth sensor 3, a fixed-star 
15 sensor 5, and an attitude angle calculation unit 8 compris- 
ing an orbit coordinate calculation unit 4, a fixed-star 
vector calculation unit 6, a yaw angle calculation unit 
7, and a command signal ON /OFF switch S. 

The orbit coordinate calculation unit 4 calculates 
20 orbital coordinates of the satellite at an arbitrary 
time instant using orbital information D transmitted 
from a ground station as a command signal. The fixed 
star vector calculation unit 6 receives data 
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from the fixed star sensor 5 and calculates therefrom a 
unit vector in the direction to a fixed star with respect 
to, the coordinates of the axis of the satellite. The yaw 
angle calculation unit 7 receives the orbital coordinates 
calculated by the orbit coordinate calculation unit 4, the 
unit vector of the fixed star under observation as 
calculated by the fixed-star vector calculation unit 6, 
and data 0(t) , 5(t) from the earth sensor 3 and calculates 
the yaw angle -+(t) of the coordinates of the satellite 
axis with respect to the orbital coordinates. 

Fig. 5 depicts a mathematical model showing the 
relation of the orbital coordinates (Xo,*0' z O) to the 
inertial space coordinates (Xi.,Yj r Zi). It is assumed that 
the ascending node right ascension, the orbital tilting, 
the period, and the ascending node separation angle at a 
time instant tg are given as the orbital information D 
from the ground command by q, i, J and Cq, respectively. 

In this case, the orbital coordinates are given 
by the following equation: 
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/ -sin f cosfl 
-cos i sinSJ cos f 



-sin i sinfl 



-sin f cosft 

+cos i cos 2 cos f 



-sin f sinfi 

4cos i oosQ cos f 



sin f cosfi 



cos f cosf; 
-cos i cosfl cos f 



cos f sin i 



-cos i 



sinn sin i 



l 0 c z ] [x x . y x , z z y 



T: transposition (2) 



f Q + 2it(t-t 0 )/T 



(3) 
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The orbital coordinate calculation unit 4 
performs its operation according to equations -2 and 3 to 
derive the orbital coordinates. 

Describing the operation of the fixed star 

vector calculation unit 6, the fixed star sensor 4 is 

mounted on the satellite with its optical axis Y s being in 

coincidence with the -Y B axis of the satellite, as shown 

in Fig. 2. A fixed star detect ion level of the sensor 4 

is set at a stellar magnitude of a star which is higher 

than third stellar magnitude, and the field of sight of 

the sensor 4 is selected as on the order of 4° x 4°. With 

these conditions, the sensor 4 mounted on a stationary 

satellite directed to the earth A can detect only Polaris 

2. Assuming that the coordinates of the fixed star 

observed by the fixed star sensor 4 are (Z,X), the unit 

vector S(f) in the direction to Polaris 2 with respect to 

Xb/ *b and Z B °f the satellite is given by the following 
equation: 

f >i ■ . ^. 

\S± f ^ cos Z sin X | 

S(f) =jS 2 = b c s cos z cos x >' . (4) 

|S 3 ^ J w sin Z : 

where (b c s! is a coordinate transformation matrix 

determined by the relation between the fixed star sensor 

coordinates system and the coordinate axis of the 

satellite/ which matrix can be obtained preliminarily. 
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The yaw angle calculation unit 7 receives IqCj] 
obtained from equation 2, the fixed star direction unit 
vector obtained by the equation 4, and the^ roll angle tf(t) 
and the pitch angle e (t) measured by the earth sensor 3 
and calculates the yaw angle * according to the following 
equation: 



sin . 
where : 



<S2bi - 33b! - S 3 b 2 + a 2 B3)/(a]bi - a 2 b2) , 



«1 ■ - '1 cos 0(t) 



(5) 



a 2 » J 2 cos fSlt) 
a 3 - / 3 C sin 0(t) 



(6) 



°1 " ff sin Mt) - / 2 C sin frf(t)cos9(t) 



'/f 

,'C 

'2 



b 2 « /, sin tf(t)cose(t) + f 2 sin 9 (t) 
b3 « / 3 ° cos ^(t)cos ©(t) 

' COS RAS COS DIS *J 

IoCiJ sin RAS cos DIS ; , ... (7) 

^ sin DIS j 

where RAS and DIS are, respectively, the right ascension 
and declination by the fixed star (Polaris 2) under 
measurement of the fixed star sensor 4, which are obtained 
preliminarily. 

According to this embodiment, the yaw angle 
can be obtained by performing operations to obtain 
Iff , 'f * ffl by inserting the output [nCj] of the orbital 
coordinates calculation unit 4 into equation 7, obtaining 
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ai and bi (i ■ 1# 2, 3) by inserting the value of 
I f\ ' h 9 fi \i thus obtained and the values (0,6) obtained by 
the earth sensor 3 into equation 6, and s calculating the 
desired yaw angle * by inserting a^ and bi thus obtained 
5 and the components (S2/S3) of the unit vector S(f) 
calculated in the fixed star vector calculation portion 6 
into equation (5) . 

According to this embodiment, the yaw angle * of 
the satellite can be obtained precisely by using the data 

10 from the fixed star sensor 4, regardless of whether or not 
the output (0, 5) of the attitude angle calculator 8 is 
nearly equal to zero or whether or not the orbital tilting 
of the satellite is nearly equal to zero. 

Fig. 6 shows another embodiment of the present 

15 invention in which the sun is used as the fixed star and a 
sun ephemeris calculator 9 is added to the attitude 
calculator 8. Other components are the same as or similar 
to those used in the embodiment shown in Fig. 4 and are 
thus depicted by the same reference numerals as those used 

20 in Fig. 4. 

With this arrangement, an orbital coordinates 
calculator 4 calculates the orbital coordinates of the 
satellite at an arbitrary time instant as described with 
reference to Fig. 4. A fixed star vector calculator 6, a 
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sun vector calculator in this case, receives data from a 
fixed star sensor, in this case a sun sensor 11, and 
calculates a unit vector in the direction toward the sun 
with respect to the coordinate axes of the satellite. A 
sun ephemeris calculator 9 operates to calculate the 
declination and right ascension of the sun in the 
direction to the sun with respect to geocentric inertial 
coordinates at the time instant at which the sun sensor 4 
detects the sun. A yaw angle calculator 7 receives the 
unit vector in the direction to the sun given by the 
orbital coordinate calculator 4, the declination and right 
ascension of the sun given by the sun ephemeris calculator 
9, and the data given by the earth sensor 3 to calculate 
the yaw angle of the satellite with respect to the orbital 
coordinates system. 

Fig. 7 depicts a mathematical model giving the 
relation of the orbital coordinate system (X 0 , v 0' z 0> to 
the geocentric inertial coordinate system (Xi,Yi,Zi). The 
right ascension « at the ascending node, the orbital 
tilting i, the period T and the ascending node deviation 
angle f 0 / all at the time instant t 0 , are here assumed as 
being given by the orbital information contained in the 
ground command. The orbital coordinates at this time are 
given by: 



-sin f cosfl -sin f sin fl cos f sin ± 

^=os i sinn cos f *cos i cos Ocos f * 9 1 



-sin i sin SI sl n f cos n -co S L 



-sin f cosn cos f cosn ^ 

♦cos i cosu cos f -cos i cosfl f 



(8) 
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where T indicate* transposition. 

I ■ fn + 2 *(t - t 0 )/T . ... (9) 

The above operations are performed .by. the orbital 
coordinates calculator 4. 

Tne operations performed by the sun vector 
calculator 6 will now be described with the assumption 
that a two-axis sun sensor of a commonly used type is 
available on the satellite. The two-axis sun sensor 
measures an azimuth angle g and an ascending vertical 
angle b of the sun with respect to the coordinate system 
of the sun sensor on the satellite, as is well known. 
Therefore, the unit vector S(t) in the direction to the 
sun with respect to the satellite axis coordinates system 
is as given below; 



f > 



S(t) 



a l cos g cos b 

?2 " IbC S J • sin g cos b i , ... (i 0) 

s 3. . sin b ! 



where ( B C S J is a coordinate transformation matrix based on 
the relation between the sun sensor coordinate system and 
the satellite axis coordinate system, which is known. The 
sun vector calculator 6 performs the above operation to 
determine the unit vector S(t). 

The operation of the sun ephemeris calculator 9 
will now be described with reference to Fig. 8, which 
illustrates a mathematical model showing the relation 



0 i 747 i 5 



- 15 - 

between the sun and the geocentric inertial coordinate 
system. In Fig* 8, the right ascension R(t) and the 
declination D (A) of the sun are given as follows: 

sin D(t) = sina sin2 (t) ... (11) 

* 

5 cos R ( A) * cos i (t)/cos D(t) ... (12) 

where a = 0 . 409 (rad) . 

t = 0.0172 x D, 

where D is the number of days counted from the equinox 
time of the year in which the satellite is launched to the 
10 current time. 

The yaw angle calculator 7 receives the unit 
vector obtained according to equation 10 , the right 
ascension and declination of the sun obtained from 
equations 11 and 12, and the roll angle #(t) and the pitch 
15 angle e (t) measured by the earth sensor 3 to calculate the 
yaw angle according to the following equations: 

sin * = <S2b, - 33b! - S 3 b2 + a 2 b 3 )/(a 1 b 1 - a2b2> # 

... (13) 

where: 
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ai = - | c cos 0{t) 
a 2 = f\ cos jz$(t) 

bi = sin H(t) - l\ sin jrf ( t ) cos H(t) ...(14) 
b 2 » /i sin #(t) cos H + / j sin H(t) 
N b 3 - cos jzJ(t) cos H(t) 

S D(t) N 

... (15) 



cos R(t) cos D(t) 
= [oCil sin R(t) cos D(t) 

. sin D ( t) 
\ ~ > \ 

According to this embodiment, the yaw angle v 

can be determined by inserting the output to c ll o£ the 

orbital coordinate calculator 4 and the outputs R(t) and 

D(t) of the sun ephemeris calculator 9 into equation 15 to 

obtain [ jf c , S\ , S\ ] , inserting [Jf C r£ c 1 and the angles 

$ and H measured by the earth sensor 3 into the equation 

14 to obtain ai, bi (i = l r 2, 3) , and inserting a^, bi 

and the component (S2FS3) of the unit vector S(t) 

calculated by the Sun vector calculator 6 into equation 

13. 

Fig. 9 shows another embodiment of the present 
invention. In Fig. 9, a register 10 stores positional 
information Sg of the satellite which is contained in the 
command signal transmitted from the ground station, and an 
initial correction calculator 12 receives the positional 
information Bg stored in the register 10 and data from the 
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earth sensor 4 to obtain an initial value of a correction 
to be added to the data from the sun sensor 4. That is, 
the initial correction calculator 12 reads out the content 
8 0 of the register 10, which gives the angular relation of 
5 the position of the satellite in a plane 1 (Pig. 3) with 
respect to the sun at the time tg and the roll angle 0(0) 
measured by the earth sensor 3 at the same time to and 
operates to obtain the initial correction ag according to 
the following equation : 

10 a 0 = (zi(O')sin 3o- . - • (16) 

A correction amount calculator 16 receives the 
data from the earth sensor 3, the initial correction ao 
from the initial correction calculator 12, and the 
angular relation between the satellite and the sun given 

15 by the orbital position calculator. 4, and from them 
calculates the amount of correction to be added 
sequentially to the data from the sun sensor 4. 

The orbital position calculator 4 reads a time 
signal At from a timer 14 and the content S 0 of the 

20 register 10 to calculate the angular relation S(At) 
between the satellite and the sun according to the 
following equation : 

6 (At) = Bo + "0 At ' • • • (17) 

where At is the time from the time instant to to and o)q is 
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the angular velocity of the satellite in its orbit (which 
is 0.25°/min when the satellite is geostationary) . 

Thus, the correction calculator 16 receives ao 
from the initial correction calculator 12 and S(At) from 
5 the orbital position calculator 4 to determine a 
correction 10(At) according to the following equations: 

10{lt) = a 1 - a 0 ... (18) 

ai = *t) sin Silt). ... (19) 

The yaw angle deviation calculator 17 receives 
10 the correction amount 10(lt) obtained according to 
equation 18, S(lt) obtained according to equation 17, and 
a variation £r ( (at) of the sun incident angle n with n (t(j} 
at time tg being zero, which is measured by the sun sensor 
4, to calculate a yaw angle deviation 1C (It) according to 
15 the following equation: 

10 (At) = [AC (It) + A0(lt) ]/cos 6(lt) . ... (20) 
As mentioned hereinbefore, according to the 
present invention, the yaw . angle of the satellite or- the 
deviation thereof can be precisely determined regardless 
20 of the roll angle $zJ, the pitch angle H of the satellite, 
the orbital tilting i thereof, and/or the angular relation 
8 between the satellite and the sun. 

Designs of components in each of the 
aforementioned embodiments may be easily made by those 
25 skilled in the art.- 
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CLAIMS 

1. Attitude angle calculation apparatus to be 
mounted on a satellite for calculating an attitude angle 
of the satellite, characterised by an earth sensor (3) 
, for detecting deviation of an angle (6(t)) defined 
5 between a predetermined axis (Zg) of the satellite and 
a straight line (PO) connecting the satellite and the 
geocentre of the earth (A); a fixed-star sensor (5) 
for detecting a physical intermatrix of a fixed star (2); 
and yaw angle determination means (4,6,7) responsive to 

10 orbital position information (D) of the satellite sent 

from a ground station and to outputs of the earth sensor 
and the fixed-star sensor to determine the yaw angle 
<*(t)) of the satellite independently of roll angle 
(<Kt)) and/or pitch angle (8(t)) thereof. 

15 2 ' Apparatus as claimed in claim 1, characterised 

in that the fixed-star sensor (5) provides coordinates 
of the fixed star in a field of sight thereof; and in 
that the yaw angle determination means comprises orbital 
coordinate calculating means (4) responsive to the orbital 

20 position information to calculate orbital coordinates 
of the satellite, fixed-star vector calculating means 
(6) responsive to an output of the fixed-star sensor 
for calculating a unit vector in a direction to the 
fixed star with respect to coordinates of the satellite 

25 axis, and yaw angle calculating means (7) responsive 

to the outputs of the earth sensor, the orbital coordin- 
ate calculating means and the fixed-star vector 
calculating means to calculate the yaw angle of the 
coordinates of the satellite axis with respect to the 

30 orbital coordinates. 
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3. Apparatus as claimed in claim -1, characterised 
in that the fixed star is the sun; and in that the yaw 
angle determination means comprises orbital coordinate 
calculating means (4) responsive to an output of the 

5 earth sensor (3) for calculating orbital coordinates 

of the satellite, a sun sensor (11) for determining azi- 
muth angle and angle of incidence of the sun with 
respect to coordinates of the sun sensor, sun vector 
determining means (6) responsive to an output of the sun 

10 sensor for calculating a unit vector in a direction to 
the sun with respect to coordinates of the satellite 
axis, sun ephemeris calculating means (9) for calculating 
the right ascension and the declination of the sun with 
respect to geocentric inertial coordinates, and a yaw 

15 angle calculator (7) responsive to outputs of the earth 
sensor, the orbital coordinate calculating means, the 
sun vector calculating means, and the sun ephemeris 
calculating means for calculating the yaw angle of the 
satellite axis coordinates with respect to the orbital 

20 coordinates. 

4. Apparatus as claimed in claim 1, characterised 
in that the fixed star is the sun; and in that the yaw 
angle determination means comprises a register (10) 
for storing positional information of the satellite and 

25 its orbit received from the ground station as a command 
signal, initial correction calculating means (12) 
responsive to the position information stored in the 
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register and an output of the earth sensor at a time 
instant at which the command signal is transmitted for 
calculating an initial value of correction, a timer (14), 
orbital position calculating means (4) responsive to the 
5 satellite position stored in the register and an output 
, of the timer for calculating an angular relation between 
the satellite and the sun, a sun sensor (5) for measuring 
a relative variation of an incident angle of the sun, . 
correction amount calculating means (16) responsive to 

10 outputs of the earth sensor, the orbital position 

calculating means and the initial correction calculating 
means for calculating a correction amount sequentially, 
and yaw angle deviation calculating means (17) responsive 
to an output of the sun sensor and the correction amount 

15 calculating means for calculating deviation of yaw 

attitude angle of the coordinates of the satellite axis 
with respect to the orbital coordinates. 
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FIG. 3 
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